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Abstract 

The impact of pionic correlations and meson-exchange currents in determining the 
(vector) response functions for electroweak quasielastic lepton scattering from nuclei 
is discussed. The approach taken builds on previous work where the Fermi gas model 
is used to maintain consistency in treating forces and currents (gauge invariance) 
and to provide a Lorentz covariant framework. Results obtained in first-order 
perturbation theory are compared with infinite-order summation schemes (HF and 
RPA) and found to provide quite successful approximations for the quasielastic 
response functions. The role of pionic correlations in hardening the responses Rl 
and Rt is investigated in some detail, including studies of the relative importance 
of central and tensor pieces of the force and of exchange and self-energy diagrams; 
in addition, their role in significantly modifying the longitudinal parity-violating 
response R^y is explored. The MEG are shown to provide a small, but non- 
negligible, contribution in determining the vector responses. 
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1. Introduction 

The role of the pion in nuclear structure has been explored for some time now ^) 
and clear signatures of its presence as a carrier of the force acting between nucle- 
ons via the one-pion-exchange potential (OPEP) have been found, for example, in 
studies of NN-scattering ^) or of the deuteron (see ref. ^) for a review in this area) 
and as a carrier of an electromagnetic current in np radiative capture ^) or in the 
electrodisintegration of the deuteron ^'^). The properties of the 3-body nuclei "^He 
and "^H also appear to be represented quite well by hadronic descriptions which in- 
clude pionic degrees of freedom in forces and currents, except perhaps at the highest 
momentum transfers in electron scattering where extremely short distance scales are 
being probed. Importantly, this success requires that the theory take into account 
2-body meson-exchange currents (MEC) and 3-body forces (see, e.g., refs. ^'^)), 
both having parts that arise from pion exchange. Such observations regarding 2- 
and 3-body nuclei support the description of nuclei in general within the framework 
of a hadronic quantum field theory (the nucleus viewed as a composite system of 
baryons — especially nucleons — and mesons — especially pions). 

However, when one comes to consider heavier nuclei the situation then becomes 
less clear because of the complexity of the nuclear many-body problem and the need 
to deal with less developed nuclear wave functions than are available for the few- 
body nuclei. To mention just a few examples where pionic effects have been sought, 
on the one hand the enhancement of the dipole photoabsorption sum rule over the 
canonical Thomas-Reiche-Kuhn value has been interpreted with some success in 
terms of a pionic contribution stemming mostly from the tensor force in second 
order ^). On the other hand, evidence for the long-sought-after softening of the 
quasielastic peak in the spin-longitudinal a q isovector nuclear response, a precursor 
of pion condensation at nuclear density, has never been convincingly found ^°'^^). 
Furthermore, while MEC effects in electron scattering from A> 3 nuclei are expected 
(see, e.g., ref. ^^)), they have been difficult to isolate from uncertainties in nuclear 
structure. 

In contemporary descriptions of many-body nuclear structure using hadronic 
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degrees of freedom different starting points have been assumed. For example, in 
the relativistic mean-field approach of Walecka and collaborators (so-called 
relativistic quantum hadrodynamics, QHD) in lowest-order the pionic degrees of 
freedom do not enter; in that case, symmetry considerations (translational invari- 
ance) prevent the pion from contributing to the energy of the system. On the other 
hand, variational approaches (see refs. ^'®)) suggest an important role for the pion: 
for example, in the work reported in ref. it is found that over half of the nuclear 
matter binding energy may be ascribed to the pion. 

As the above issues are not yet fully resolved, we are left with some uncertainty 
about the size of the role played by pionic degrees of freedom in many-body nuclear 
structure. When attempting to elucidate this issue, it should be kept in mind that 
these degrees of freedom may be of rather different importance in describing different 
nuclear observables. In particular, their role in accounting for the binding energy 
of nuclei may be rather different from that involved in describing the electroweak 
quasielastic nuclear response functions which are the main focus of the present work. 
In the former, (particle-particle) interactions between nucleons in the Fermi sea are 
presumed to provide most of the binding energy, whereas in the latter particle-hole 
excitations involving particles typically at several hundred MeV in the continuum 
are dominant. 

In this paper, which enlarges on a previous one ^^), we do not aim at establish- 
ing a relativistic quantum field theory of nuclei based on hadronic degrees of freedom 
or to explore ground-state properties using sophisticated nuclear wave functions. 
Rather we treat pionic degrees of freedom as perturbations, building on a simple, 
tractable, covariant nuclear model for quasielastic response functions, namely, the 
relativistic Fermi gas (RFG) model ^^). The free RFG model appears to describe 
nuclear excitations at high three-momentum transfer q and high energy-transfer 
u) reasonably well for kinematics near the quasielastic peak (QEP), that is, where 
uj ^ I (5^ I /2m AT, with = uP' — . To the extent that pionic correlations and 
pionic MEG provide only moderate corrections to the free RFG responses, one may 
hope that a perturbative treatment is meaningful. Accordingly, we start with all 
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Feynman diagrams carrying one pionic line and later extend our analysis to account 
for perturbative diagrams up to infinite order to test their importance. Issues re- 
lating to Lorentz covariance and electroweak gauge invariancc can be explored in 
this model where a high level of consistency in treating forces and currents can be 
maintained. Some of the issues relating to gauge invariance in this pionic model 
are addressed in the present work. In subsequent work it will be important to ex- 
plore the roles played by mesons heavier than the pion within the same framework, 
however, the present scope has been restricted to a study of the pionic effects in 
quasielastic electroweak nuclear response functions. Given the long-ranged charac- 
ter of pion exchange we believe that such effects provide a natural starting point 
for a more ambitious study. 

In the present work we study the pionic correlation efi'ects embodied in the 
so-called self-energy and exchange diagrams, whereas, for the MEG contributions 
we shall mostly rely on past work ^^), paying special attention to the problem of 
fulfilling the continuity equation in achieving a consistent treatment of currents and 
forces. In accord with those previous studies, we find sizable pionic contributions 
to the electromagnetic longitudinal (spin scalar, cr = 0) and transverse (spin vector, 
cr = 1) nuclear responses. In both cases it is found that the correlation effects 
produce a hardening of the responses, that is, a shift of the strength to higher w. 

In particular, as discussed in detail in the following sections, from our past and 
present analyses the following points emerge: 

a) In isospace the contribution of the self-energy diagram to the charge response 
is almost equally split between isoscalar (r = 0) and isovector (r = 1) com- 
ponents. (The latter, on the other hand, is of course overwhelming in the 
transverse response, due to the dominance of the isovector magnetic moment.) 
In contrast, the r = part of the exchange diagram is three times as large 
as the r = 1 one in the charge response and this imbalance, which becomes 
even stronger in higher orders of perturbation theory, is further strengthened 
by the difference between the isoscalar and isovector form factors (see later). 
The isoscalar dominance of the pionic exchange correlations has dramatic con- 
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sequences for the weak neutral current fongitudinal response function, as we 
shall see in sect. 4. 

b) While the tensor component of the OPEP never contributes to the self-energy 
in a translationally invariant system, the exchange diagram gets a tensor contri- 
bution, but only in the transverse channel and mostly via the backward-going 
graphs. This implies a different role for the pionic force in the two electromag- 
netic responses, a finding which should be tested against experiment. 

c) When we extend our analysis from first to infinite order of perturbation theory, 
thus generating the Hartree-Fock (HF) approximation with the self-energy dia- 
gram and the random phase approximation (RPA) with the exchange diagram, 
our results do not change substantially. It therefore appears that, although the 
pionic interaction is strong, nevertheless for quasielastic kinematics its effects 
are reasonably small at high-g, thus rendering perturbation theory quite ac- 
curate already at the lowest order, at least for the classes of diagrams studied 
here. 

d) Finally, the contribution of the central component of the pionic interaction to 
the exchange diagram stems largely from the 5-force and not from the finite- 
range one. In keeping with the usual approach taken in studies of pionic effects, 
we include a ttNN vertex function F^r, whose scale is set by a mass parameter 
Att, to smear out the 5-piece of OPEP. Consequently, one of the goals in the 
present work is to explore the sensitivity of the resulting nuclear quasielastic re- 
sponses to the choice of the phenomeno logical parameter A^r. The question then 
arises: To what extent is the continuity equation modified by the presence of 
F^? We address this issue later in the present work and here simply recall that 
for pointlike nucleons the continuity equation is indeed obeyed by the OPEP 
and by the related pionic MEC ^^), but is not affected by the 5-interaction. 
This conclusion must be modified, however, when the vertex function F^^ is in- 
troduced, although, as discussed below, additional MEC can be introduced 

in a way such that the continuity equation is still satisfied. 
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In contrast to the exchange diagram, the self-energy term gets a contribution 
only from momentum-dependent forces, and therefore an unmodified pionic S- 
interaction does not contribute in this channel; in fact, the contribution coming 
from the 5-function in OPEP, when modified by the ttNN vertex form factor, 
contains an effective momentum-dependence and so is nonzero. 
Regarding the two parameters that characterize our approach, namely the 
Fermi momentum kp and the mass parameter Att, we note that the value of the 
first is essentially set by the nuclear density. Ultimately, given extensive evaluation 
of all nuclear model dependences (correlation effects, MEC effects, etc.), electron 
scattering measurements in the quasielastic region should serve to determine kp- 
Likewise, some information on the range of acceptable values for the phenomenolog- 
ical parameter can in principle be obtained from comparison with experiment. 
Clearly the value deduced for this parameter is model-dependent: it is meant to 
absorb effects from hadronic physics other than the explicit pionic degrees of free- 
dom (to the extent that this is even possible through a single vertex function F^). 
In an extended model with other active mesonic degrees of freedom presumably 
the physics embodied in this way can be quite different. Even given the limited 
scope of the present pionic model, it is interesting to explore the sensitivity in our 
results to the actual choice of A^r and we do so in sect. 5. The unexpected (and 
perhaps important) finding in this connection is the following: the diagrams which 
contribute most to the nuclear responses are those that are least affected by F^. 

With respect to the relativistic aspect of the present approach, one limitation 
of our treatment arises from the static character of our pionic interaction, while 
another arises from the not-fuUy-relativistic character of the fermion propagators 
that we use. However we have been able to achieve an almost exact relativistic 
treatment of the fermion kinematics. As discussed in ref. ^^), this is obtained 
by expressing the electromagnetic longitudinal and transverse responses and, as 
well, their weak neutral current analogs in terms of a single relativistic scaling 
variable ^^). Moreover we use an approximation to the relativistic electromagnetic 
and weak neutral current vertices, which has proven to be quite accurate. As we 
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shall see, this result is of relevance because the pionic effects in these observables 
are felt up to quite large momentum transfers (~ 1 GeV/c) where relativity cannot 
be ignored. 

The items elaborated upon in the present paper are treated in the following 
order: in sect. 2 we deal with the self-energy diagram both in first-order pertur- 
bation theory and in the HF scheme. In sect. 3 we consider the exchange diagram, 
again both in first- and infinite-order perturbation theory, where the latter is ex- 
plored within the framework of continued fractions. In first-order, both forward- 
and backward-amplitudes are included in the analysis. In sect. 4 we calculate the 
electromagnetic and the weak neutral current responses, focusing in particular on 
the large enhancement found for the weak neutral current longitudinal response 
when correlations are included, compared with the results found with the free RFG 
model where a delicate cancellation occurs. In the same section we also address the 
question of the evolution with q of both the correlation and MEG contributions. 
In sect. 5 we study the At^- and A;i?-dependences of our results. In the concluding 
section we summarize our results as they stand at present in treating pionic effects 
and their impact on quasielastic nuclear response functions. 



- 8- 

2. Self— Energy Contributions 

In this section we explore the contribution of the self-energy diagrams to the longitu- 
dinal and transverse electromagnetic responses of an infinite, homogeneous nuclear 
system with an equal number of protons and neutrons (Z = N = A/2). The corre- 
sponding (vector) weak neutral current responses are considered in sect. 4, whereas 
the axial-vector weak neutral current response is studied in a companion paper ^^). 
We shall confine ourselves to a consideration of the one-particle-one-hole (Ip-lh) 
sector of the excitation spectrum of the RFG and, as already stated in sect. 1, in 
our scheme only the pion mediates the interaction between the nucleons. 

The nuclear responses may be expressed via the polarization propagator. In 
first-order perturbation theory the pion yields two self-energy contributions to the 
latter, by dressing either the particle or the hole propagation in the associated 
Goldstone diagrams. These are displayed in fig. 1 and the corresponding analytic 
expressions for both the longitudinal and transverse channels have been derived in 
ref. 1^): 

Ai;-,_,fe.)=C.,.^lim|- 

(27r)3 ^ ^^l '^^^ ^(p_^)2_^^2 ttvi (^q k - p)^ + ml ) 

(2.1a) 
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Here the contribution associated with the convective nucleon current has been ne- 
glected, since it represents only a tiny fraction (a few percent, at most) of the total 
transverse response for momentum transfers exceeding about 300 MeV/c. Note 
that in (2.1) the derivative and the limit with respect to the variable a define the 
prescription one should follow in order to deal with the double pole characterizing 
the first-order self-energy contribution to the nuclear responses. 

In the above, along with the Fermi momentum kF, the three-momentum trans- 
fer q, the energy transfer uj, and the spacelike four-momentum transfer = 
itp' — < 0, with q — |g|, the dimensionless variables of ref. have also been 
employed: 



K, = 



q 



2m N 



X = 



2m jv 



riF = 



T = 



rriN 



4:m%' 



rriN 



and = 



mN 
(2.4) 



Moreover, //^ = //p — and = fip + fin are the isovector and the isoscalar 
magnetic moments of the nucleon, respectively. To incorporate some aspects of the 
ttNN vertex, we employ the monopole form 

Al - ml 



(2.5) 



A2 +p2 ' 

which will be discussed at some length later on, especially in connection with the 
problem of choosing an appropriate value for the phenomenological constant A^r. 



Notice that in the energy-conserving (5-function occurring in our expressions 
for the self-energy contributions to Rl and Rt we have replaced the non-relativistic 
term q'^/2mN with the relativistic one \Q^\/2mN- This can be shown to be al- 
most exactly equivalent to employing relativistic kinematics and leads to the scaling 
variable of the relativistic Fermi gas model which enters in (2.1): 



V'r = — 

Vf 



A(A + 1) 

K, 

K 



(2.6a) 



(the last terms on the right-hand side of (2.1b) correspond to the Pauli-blocked 
region, where no scaling occurs). 

Had we kept the term j^m^ in the energy-conserving 5-function, then we 
would have obtained the same results but with the non-relativistic scaling variable 

V'nr = —(--«) (2.66) 

rjF / 

rather than i/jr- Actually, an additional factor 1+2 A also appears in the "relativized" 
responses: it expresses the Jacobian of the transformation from the variable ipnr to 
V'r and has been included in the multiplicative form factors (see below). 

Another place where relativity is carefully accounted for in our approach is 
in the electromagnetic vertices. Indeed, the electromagnetic form factors /l,t(Q^) 
presently employed are given by the following expressions: 

/E(g^) = (l + 2A)|^[G^nr)]' + r[G^*(T)]2^(l-V',2)} (2.7a) 

/|(g2) = {l + 2X)j-^[GZ'{r)f, (2.76) 

with 

G^*(r) = C-+mt)GE,{r) + (^ - mOG^^ (r) (2.8a) 

GZ'ir) ^ C-+mt)GM,{r) + {-^ - mt)GMAr), (2-86) 

Gep „ and Gmp „ being the Sachs form factors of the proton and neutron and 
mt labeling the isospin projection. The results in (2.1) are actually obtained by 



adding the proton response {mt = 1/2) multiplied by Z to the neutron response 
{rrit = —1/2) multiplied by N. Formulae (2.7) and (2.8) are deduced in ref. and 
there shown to provide an accurate representation of the RFG responses. We recall 
that the need for a fuUy-relativistic treatment of the 7NN vertex arises from the 
existence in the problem of the non-relativistic reduction of the nuclear responses 
of a second scale involving n and the squares of the nucleon magnetic moments ^^), 
beyond the usual one set by k alone. 

A final comment on the expressions (2.1b), already discussed in ref. ^^), relates 
to the domain in the (g,a;)-plane where AWj^^J"^ ^ and AR^^J'^ ^ are nonzero. It is 
the same as that of the RFG, although the self-energy contributions do not vanish 
on the boundaries of the response region and are discontinuous across the boundary 
dividing the Pauli-blocked region from the non-blocked one (see the case at q = 300 
MeV/c below in fig. 2). This is connected with the fact that the derivative of the 
RFG response does not vanish on the boundaries themselves and is discontinuous 
as well when one crosses into the Pauli-blocked region. 

As already mentioned, in the present approach we treat the pionic lines of fig. 1 
at the static level and hence they correspond to the well-known one-pion-exchange 
potential (OPEP): 



In (2.9) one easily recognizes the central contact and momentum-dependent pieces 
of the interaction in addition to the tensor contribution. By performing spin traces 
(see Appendix A), it is then an easy matter to show that the tensor force S12 
never contributes to the self-energy diagrams of fig. 1, i.e., in either longitudinal or 
transverse channels. 

It is equally easy to show, by the same token, that the longitudinal isoscalar 
(r = 0) and the isovector (r = 1) self-energy contributions are identical, but 
for the difference between the isoscalar and the isovector form factors and 



V{q) 




m: 



2 



(2.9) 



fl . This difference arises on the one hand because the neutron charge form 
factor is nonzero (leading to a small effect). On the other hand, the expressions 
(2.7a) contain terms which involve the proton and neutron magnetic form factors 
and these conspire to yield a much larger isovector than isoscalar contribution (see 
below). Since these GM-effects are multiplied by r]p ^ 0.04-0.08, they tend to 
be suppressed; however, since they are multiplied by r, they provide increasingly 
important contributions as the momentum transfer increases. The usual approach of 
including these (relativistic) corrections via the Darwin-Foldy term already breaks 
down for g < 1 GeV/c, as discussed in ref.^^). 

From (2.1) it clearly appears that the self-energy contribution actually results 
from the cancellation of two pieces associated, respectively, with the dressing of the 
hole and particle propagation: the larger term is the latter (at least for not too 
small and q not too large, see sect. 4 and fig. 2). Therefore, Ai?^°'J^ depletes the 
nuclear response at low-a;, while enhancing it at high-a;, as it should in order to 
obey the sum rule requirement ^^). Worth noticing is that the contribution of the 
contact (5) term of the OPEP cancels in (2.1), but for the momentum-dependence 
introduced by the vertex function Ft^. When the latter is taken into account, as in 
our case, then an additional cancellation between the S- and momentum-dependent 
contributions occurs. What remains is then a rather modest contribution (negative 
at small and positive at large uj, again for not too small and not too large q), 
that rapidly decreases as q increases. This is illustrated in fig. 2 for q = 300, 500 
and 1000 MeV/c in the longitudinal channel. For sake of illustration here and in 
the following we take kp = 225 MeV/c (which roughly corresponds to the case of 
^^C) and Att = 1300 MeV, unless otherwise specified. 

A similar situation also holds in the transverse channel, but for the well-known 
factor Zfip + Nji^ = \A{ji'^ -\- jj^) = |^(//| + A*^)) which greatly suppresses the 
isoscalar contribution compared to the isovector contribution because of the fact 
that ins/l^v)^ = 0.035. 

An issue we wish to address in closing this section relates to the self-energy 
diagrams of higher-order in perturbation theory. Do they play an important role? 



To answer this question at least partially, we may examine the longitudinal electro- 
magnetic response in the Hartree-Fock approximation (actually here only the Fock 
approximation, since the Hartrcc contribution vanishes in nuclear matter, although 
we shall continue to use the label HF) and compare with the results displayed in 
fig. 2 to which we add the free response. The HF longitudinal response reads 



is the pion self-energy. 

In fig. 3 the longitudinal HF response is displayed for q =300, 500 and 1000 
MeV/c. Note that this response appears in a region of the (q',a;)-plane somewhat 
different from that where the RFG responses are nonzero, notwithstanding that each 
order of perturbation theory in the self-energy is actually nonzero precisely where 
the RFG is nonzero. The effect, however, is so small (at most a few MeV) that it 
can only be perceived at \ow-q, as seen in fig. 3. Remarkably, from this figure one 
realizes that the first-order response practically coincides with the HF result: the 
many-body framework apparently quells the strong interaction carried by the pion 
to the point of rendering quite accurate a first-order perturbative treatment of the 
OPEP, at least as far as the mean-field contributions are concerned. Accordingly, 
the impact of the pion self-energy contribution on the total charge response is 
rather mild and one observes some hardening for q 300 MeV/c, which, however, 
disappears as q increases (see sect. 6). 




(2.10) 



where 




(2.11) 
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3. Exchange contributions 

In the Ip-lh sector of the nuclear response, in addition to the self-energy diagrams 
dealt with in the previous section, two further classes of Goldstone diagrams yield 
contributions in first-order perturbation theory: they are displayed in fig. 4 and 
are commonly referred to as exchange diagrams. In particular, in fig. 4a the Ip-lh 
state propagates forward in time (as in the so-called Tamm-Dancoff approximation 
(TDA)), whereas in fig. 4b a backward-going Ip-lh state occurs (as when one incor- 
porates ground-state correlations in the RPA). The forward- and backward-going 
first-order perturbation theory contributions will be denoted F and B, respectively. 
The expressions for the associated contributions to the electromagnetic charge re- 
sponse read: 

2 P. 



U3 



mi 



[ + q\- kF)e{kF - ki)e{\k2 + q\- kF)e{kF - k2) 
oj-\Q^\/2mN -q-k2/mN ml + [ki - k2r 



1 rl 



x l9{r]F - n) I dki j c?A;2 [A/l(/ci, A;2; A;2, V'r) - -A/K/ci, A;2; t/^r)] 

I Jk< J\2K/r]F-l\ 



1 rk> 



+ d{2K-r]F) / dki / dk2[J^L{ki,k2;k2,ipr) - J^L{ki,k2;-k2,ipr)] 
Jk< Jo 



(3.16) 
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dk2 [Miiki, k2; k2, -^'r) - A/l(A;i, /c2; t, -ijjr)] 



(3.26) 



respectively. 

In the transverse channel we have instead 



+q\- kF)d{kF - k^)e{\k2 + q\- kF)e{kF - k2) 

(h ^ ^(^- IQ^I/2mjv-g-fci/mjv) [9 ■ (fci - ^2)]^ 



Vf ruN 



xle{r]F-n)[ dki [ dk2[NT{ki,k2\k2,'^r) -■hfT{ki,k2\t,'^r)\ 

y Jk< J\2K./r]p-l\ 

+ d{2K-r]F)MWhq<)-nTik^,q<)] \ (3.36) 



and 



X J (^^(1^1 + kF)e{kF - ki)e{kF - \k2 + q\)e{k2 - kp) 

a;- |(52|/2mjv-g-K2/mjv ml + [ki - k2)^ 
x<e{r]F-K) dki dk2[- MT{-kuk2;k2,-'4)r) +J^T{-kuk2\t,-'4)r)\ 

Jkf J\2K/'nF-l\ 

-e{2n-i^F)^r[nT{k>,q<)-nT{k<,q^)\\ . (3.46) 
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In (3.1-4) the following quantities have been introduced: 
t = 



N 



k< = 



^2 



\lpr\i Ipr > — 2K/r]F 

^k/tjf — 1, K < rjF 
1, K> rjF 



= Vl + {'^l^hpf + 4:Klpr/VF = ^1 - {^k/tIfY + AKipr/VF 



9< = 



(2K/'r]F, K<riF 
\ 2, K>'qF 



whereas the functions Nl,t and Tix are defined in Appendix B. Note that while 
we have been able to express the exchange contribution to the transverse response 
analytically, at least when k, > r]F, this turned out not to be possible in the case of 
the longitudinal channel. 

A few remarks should be made concerning the above expressions. First, notice 
that all of the exchange contributions are different from zero in the response region 
of the RFG, vanishing, however, on its boundaries (at variance with the self-energy 
case). Second, we note that the exchange contributions can be split, as in the self- 
energy case, into isoscalar and isovector components. By carrying out the pertinent 
traces over the Pauli matrices (see Appendix A), in this case it is found that in 
the r = channel the pion exchange correlations are three times stronger than in 
T = 1 case, being moreover of opposite sign in the two isospin channels. Finally, 
performing a similar analysis in spin space (see Appendix A), one finds that the 
tensor interaction 5'i2 does not contribute to the charge response, but only to the 
spin response. 

The results of our calculation of the exchange contribution to the nuclear re- 
sponses are displayed in fig. 5, where both the forward- and backward-going ex- 
change contributions to the charge response are displayed for q =300, 500 and 1000 



MeV/c, and in fig. 6, where the same is done for the spin response. A feature 
common to all of the results we present here is immediately apparent: strength is 
removed by the picnic exchange correlations from the low-energy side of the nuclear 
responses and added to the high-energy side. This hardening effect goes in par- 
allel with the one induced by the self-energy correlations discussed in sect. 2, which 
however was somewhat milder. As a consequence both charge and spin responses 
turn out to be hardened with respect to the corresponding free RFG responses. 

The hardening of the charge response arises from the strong repulsive character 
of the exchange isoscalar pionic correlations, which are overwhelmingly carried by 
the (5-component of the OPEP (the attractive isovector correlations for zero-range 
contributions are three times weaker, as mentioned above). The finite-range central 
interaction of the OPEP, while attractive is not very effective and only modestly 
reduces the impact of the 5-force on the charge response. Also worth mentioning 
is that the role of ^-contributions is clearly apparent in the F-term and becomes 
dominant in the S-term. 

In the spin channel the isoscalar correlations are dramatically suppressed rela- 
tive to the isovector correlations by the factor = 0.035. In contrast to the 
situation for the charge case, the isovector central correlations arc now of a repul- 
sive character (see Appendix A for the derivation of the related spinology), leading 
again to a hardening of the spin response due to the action of the 5-force. The lat- 
ter, however, is counteracted here not only by the attractive momentum-dependent 
central interaction, but by the tensor force as well. Accordingly, the amount of 
hardening of the spin response is somewhat moderated with respect to the charge 
case. 

Interestingly, in the spin channel the (5-dominance is restricted only to the 
F-diagram, since the largest contribution to the -B-diagram is in fact provided 
by the tensor force — and the latter is known to be quite efficient in providing 
ground-state correlations. This explains why the contribution from the S-diagram 
survives in the spin channel up to larger momenta than in the charge channel. With 
regard to the momentum-dependence of the exchange contributions to the nuclear 



responses, one should also notice that the S-term falls faster than the F-term as 
the momentum transfer increases because, as seen from (3.1b), (3.2b), (3.3b) and 
(3.4b), the integrals involved in the former span a higher momentum range than 
those in the latter. As a consequence, the backward-going responses are cut-off 
more rapidly by F^r than are the forward-going ones. 

In summary, from the above analysis the following observations can be made: 

a) the exchange diagrams of fig. 4 are the ones that contribute most strongly to 
the nuclear responses as first-order perturbations in the Ip-lh sector of the RFG; 

b) their magnitude and sign is to a large extent set by the 5-component of the 
OPEP. 

This last finding will be addressed further in sect. 5 in connection with the discussion 
of the global gauge invariance of our theory. 

A further aspect of our results worth noting is that the sum-rule requirement 
is indeed fulfilled by the F-contribution, whereas it is violated by the S-term. ^^) 
Indeed, the former is characterized by a node in the response function (when ex- 
pressed as a function of a;) located at the peak position of the RFG response. Here 
the very tiny discrepancies seen can be ascribed to the present approximate, al- 
though very good, treatment of the relativistic kinematics. In contrast, the latter 
never changes sign, always remaining negative. In this connection one sees that the 
hardening of the nuclear responses previously discussed clearly follows from the way 
the sum rule is fulfilled by the F-diagram. In addition a significant contribution to 
the hardening also comes from the 5-diagram. In fact, because the (negative) B- 
term obeys the energy-weighted sum rule, it then reaches its minimum somewhat 
before the peak of the free response. 

To conclude this section we now explore the role played by exchange diagrams of 
order higher than one in the nuclear response. We do so for the charge longitudinal 
channel in the TDA, which amounts to accounting for all terms in the series shown 
in fig. 7. To achieve this, we employ the method of continued fractions ^^): in this 
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framework the polarization propagator becomes 

uIKk, A) = , (3.6) 

^ i + nW>,A)/no(«,A) + ... ' ^ ^ 

where is the longitudinal first-order exchange contribution in the isospin chan- 
nel T and n° the RFG propagator. The imaginary part of proportional to 
(3.1), whereas the real part is given in Appendix B. 

We limit our attention to the first one of the infinite set of continued fractions 
embodied in (3.6), recalling that in the limit of a zero-range force the first iteration 
of (3.6) actually corresponds to the exact TDA solution. As previously discussed, 
since the 5-piece of the interaction plays the largest role in the charge response, we 
believe that the results displayed in fig. 8 are quite representative of the true TDA 
solution (of course, one should also take into account that in coordinate space 
actually smears out the 5-f unction). 

In order to characterize the results obtained here, in fig. 8 two cases are shown: 
i) the longitudinal TDA exchange response with a pure (^-interaction together with 
and ii) adding to this the momentum-dependent forces as well. The latter, which 
corresponds to a weaker interaction, shows that the first-order result seems to be 
quite accurate even at q =500 MeV/c. However, the degree of accuracy appears to 
be lower than in the corresponding self-energy case, in accord with the finding that 
the action of the pion is stronger in the exchange channel than in the self-energy 
channel. 



- 20- 

4. Electromagnetic and weak neutral current responses 

In the previous two sections we have set up all the ingredients needed to calculate not 
only the electromagnetic responses, but also the weak neutral current ones 23-27^ 
first- (and as well in infinite-) order perturbation theory, within a framework where 
only the pion-induced correlations are taken into account. Indeed, by splitting the 
electromagnetic responses into their isospace components according to 

= i?|-(r = 0) + W^{t = 1) (4.16) 

and exploiting the conservation of the vector current (CVC) ^^), one easily ob- 
tains the weak neutral current longitudinal and transverse responses (as mentioned 
earlier, the axial-vector one is discussed in a companion paper ^^)). They read 

R'iy = aA [p^v^RTir = 0) + p'y'^Rrir = 1)] (4.2a) 
Rly = aA [p^v^RTir = 0) + P^y^RTir = 1)] , (4.26) 

where oa = — 1 and the isoscalar and isovector weak neutral current couplings, Py^ 
and Py^ , respectively, are fixed at tree-level in the standard model to the values 

= -2 sin^ Ow ~ -0.454 (4.3a) 

and 

/jW = 1-2 sin^ 9w ~ 0.546 , (4.36) 

using sin^ 6w ~ 0.227 for the Weinberg angle. 

To the contributions stemming from the correlations, those arising from the 
MEC should then be added. These have been thoroughly discussed within the con- 
text of Ip-lh excitations in ref. ^^), where their contribution to the electromagnetic 
responses has been calculated including terms of order up to in a non-relativistic 
reduction, where, in addition to spatial components, the time component of the 
MEC has also been kept in the evaluation. The ensuing violation of the continuity 
equation thus introduced has been checked in ref. ^^) and found to be quite small. 



The associated expressions in the weak neutral current sector are most easily ob- 
tained by recalling that the MEC for leading-order 7r-exchange are purely iso vector 
in nature. Accordingly one simply obtains 



(4.46) . 



(4.4a) 



In figs. 9 and 10 we now display the electromagnetic responses (longitudinal 
and transverse, respectively) at q =300, 500 and 1000 MeV/c. We observe that in 
both channels a hardening of the responses occurs, particularly at low-g, essentially 
induced by the 5-force via the exchange diagram, as discussed at length in the pre- 
vious section. Also in accord with the considerations developed there the shift in 
the maximum of the response is seen to be larger in the charge channel than in the 
spin one. As the momentum transfer increases, the amount of hardening decreases, 
since the pion-induced correlations weaken for large q. This behaviour partly re- 
flects the role of the hadronic vertex function in cutting out the high-momentum 
contributions to the response functions and partly reflects the constraints arising 
from the many-body system itself, which appear through the mismatch between the 
single-particle wave functions of the particle-hole pairs involved in the intermediate 
states of the polarization propagators considered here. 

Clearly, a unique way to assess the impact of correlations on the response 
functions does not exist. Among several different possibilities, we have chosen to 
focus on the following two quantities: i) the shift Aa;QEP of the peak position 
of the response with respect to the corresponding quantity in the RFG, which is 
displayed versus q in fig. 11 for both the longitudinal and transverse channels; ii) 
the correlated sum rule, i.e. the integral of the nuclear responses over uj, divided, 
in the longitudinal and transverse channels, by the functions 



Af 



[(1 + t)W2{t) - Wi{t) + W2{t)A] 



for L 



Cl,t{t, A;vf) 



X 



T 



for T 



2Wi{t) + W2{t)A 



(4.5) 
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respectively. As shown in ref. ^^), this optimally separates the nucleonic physics 
from the many-body problem for the RFC In (4.5) 

Wi{t)=tGI,{t) (4.6a) 

^2^^) = ^^'^^^^ + rGlir)) . (4.66) 

and 

A = ^(1 - 1^,^!^+ i(l - V^,')^^! , (4.7) 

ipr being the scaling variable (2.6a) and = ^Z^ + Vf ~ ^- recall that with 
these definitions we actually integrate the part of the responses which scales. So 
it should not be surprising that the asymptotic value, namely one, is reached for 
momenta much larger than 2k p- The normalized sum rules are reported versus q 
in figs. 12a and 12b, again for both channels. We recall that the sum rules refiect 
only the roles of the backward-going correlations and of the MEC. 

In accord with the discussion of sect. 3, from these figures one first infers that 
the correlations persist up to larger momenta in the longitudinal channel than in 
the transverse one; secondly, one sees that the MEC contribution, although not 
large, is significant, rather constant with q and acts in a way that depletes the sum 
rule. Finally inspection of the responses shows that in the longitudinal channel the 
S-terms decrease in magnitude with q more rapidly than the correlations stemming 
from the F-diagram, whereas in the transverse channel the opposite occurs — a 
clear signature of the action of the tensor force. 

Turning now to a discussion of the weak neutral current responses, we first 
consider the charge channel. To understand the nature of the problem, we display 
in fig. 13 the longitudinal nuclear response at q =300 MeV/c separated into r = 
and T = 1 channels. There we see that the action of the pion is very important in 
hardening the isoscalar channel, while rather gentle in softening the isovector one. 
Note that the effect of the longitudinal electromagnetic form factors in contributing 
to this imbalance is negligible at this low momentum transfer. When we combine the 
two isospin contributions with weighting factors ^ and Py^ according to (4.2a) 
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to get the longitudinal weak neutral current response we obtain the result shown in 
fig. 14. While the magnitude of free response is quite small, the pion-correlated one 
is dramatically modified. The origin of this effect was discussed previously in ref. 
(see also refs. 27,24^^^ where it was noted that the delicate cancellation that leads to 
the suppression of R^y in the RFG model is broken when the isoscalar and isovector 
channels are correlated differently. This heightened sensitivity to isospin correlation 
effects opens a new window on nuclear physics in the quasielastic region, possibly 
offering the unique possibility of separating the two isospin channels, something that 
is impossible to achieve with parity-conserving electron scattering alone. It can, 
however, be brought to light using parity-violating electron scattering, as discussed 
in ref. ^^), where it is shown that isospin correlations have dramatic consequences for 
the forward-angle parity-violating asymmetry over a range of momentum transfers 
persisting up to g ~ 500 MeV/c. 

In concluding this section, we return to touch on the transverse weak neutral 
current response R^y- ^^^^ factor suppresses the isoscalar contributions and 
consequently no sensitive cancellation occurs in the free RFG model to be broken 
by isospin correlations, in contrast to the case of R^y- Accordingly, the prediction 
for R^y is very simple: it is approximately half {Py^ = 0.546) of the corresponding 
electromagnetic response (of course in (4.2a,b) the electromagnetic and weak neutral 
current coupling constants are not included). These observations were critical in 
our earlier work ^^) where the insensitivity of the transverse response functions 
to isospin correlation effects allowed us to suggest backward-angle parity-violating 
quasielastic electron scattering as a tool to probe the single-nucleon form factors 
themselves. We shall return to quantify the roles played by pionic effects in this 
situation in ref. ^^). 



5. and kp dependence of the nuclear responses 

In this section we explore the significance of the roles played by the two parameters 
on which our calculation of the nuclear responses depends, namely and kp. We 
start by considering A^, which incorporates some aspects of the short-range physics 
affecting the pionic correlations. To get a feeling for the distances involved in our 
problem, we consider the Fourier transform of F^r (see (2.5)), namely 

F.(r) = (A^-m^)^. (5.1) 

The corresponding rms radius of the interaction region 

, r r -|i/2 ,/a , 

\/<7^ = [ / r^F, (r )J = j- Vl - {mJA^)^ , (5.2) 

is then easily obtained. Typical values of A^^ are of order 1-2 GeV and hence 
(m7r/A^)^ ~ 0.005-0.02, namely, only a small correction in (5.2). Accordingly the 
rms radius above is given to a good approximation by -s/G/Att: for the value A^r = 
1300 MeV used in most of this work, this yields an rms radius of 0.37 fm. 

On the one hand, as seen in fig. 15 the sensitivity of the exchange contribution 
to variations in A-^- appears to be quite mild. For sake of illustration we limit 
ourselves here and in the following to consideration only of the q = 500 MeV/c case 
and only of the longitudinal channel. The self-energy contribution, on the other 
hand, is more sensitive to changes in A^, even reversing its behaviour in the energy 
variable for A^ ^ 800 MeV, reflecting the weakening of the particle self-energy with 
respect to that for the hole (fig. 15). 

Somewhat in between the two situations discussed above is the behaviour versus 
Att of the MEG (fig. 16). We may conclude that, owing to the major contribution 
to the nuclear responses arising from the exchange term, our pionic model is indeed 
affected, but only mildly, by the short-range physics embodied in the vertex function 
F^ for A^ > 1.3 GcV. 

Glearly, all of the short-range physics cannot be represented solely by the 
monopole form factor in (5.1). An adequate treatment of the short-range physics 
(presuming we continue to adopt a strategy of sticking to a mesonic model of the 



nucleus) would require on the one hand the introduction of heavier mesons, and 
on the other the calculation of additional contributions (e.g., the ladder diagrams) 
beyond the self-energy and exchange terms considered in the present approach. 
An alternative option would of course be to abandon the hadronic description of 
the short-range physics in favour of QCD degrees of freedom, but this has not yet 
proven to be tractable in the strong-coupling regime. 

Next let us comment on how the continuity equation is affected by the vertex 
function Fj^. In this connection we recall that the continuity equation in momentum 
space reads 

i iPi +P2) ■ Jmec (Pi,P2) = (ri X [Vip2) - V{pi)] , (5.3) 

where Pi and P2 = Q — Pi are the momenta carried by the meson. For pointlike 
nucleons (5.3) is fulfilled by the OPEP (2.9) and by the longitudinal components of 
the pion-in-flight and contact currents: 



f2 

{Pl +P2) ■ Jtt (Pl,P2) = {Ti X T2) 



1 1 
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(Pl +P2) • 3 contact (Pl,P2) = (^^1 X '^2)^ 

ai ■ Pi£T2 ■ P2 ^ ^1 P2^2 ■ P2 CTl PlCf2 " Pl ^1 Pl^2 P2 



(5.46) 



+ p^ + p^ + pf + pf 

It is of significance that the 5-piece of the OPEP does not contribute to (5.3), 
since its commutator with the charge operator vanishes. Accordingly, to keep or 
to drop the 5-term in the potential is of considerable relevance for the nuclear 
responses, but is irrelevant as far as the global gauge invariance of the theory is 
concerned. It becomes of importance, however, when the vertex function is 
included. In fact, in order to satisfy the continuity equation in such a case new 
terms should be added to the pion-in-flight current, for example, the following ^^) 
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(5.5) 



{Al+pl){ml+pl) 

In contrast, no modifications are necessary but for the multipficative factor 
occurring in the contact (or pair) term, this being the one contributed to by the S 
(smeared) piece of the OPEP. The first and the third terms on the right-hand side 
of (5.5) may be viewed as describing the couphng of the virtual photon to a fictitious 
particle carrying the same quantum numbers as the pion, which thus underlies the 
microscopic description of F^r. We have numerically checked the importance of 
these additional terms and found them quite small (about 10% of the pion-in-fiight 
contribution at the peak of the latter in the charge channel and about 20% in the 
spin channel under the same conditions and so yielding less than 1% corrections 
to the total responses). This occurs unless becomes unreasonably small — say 
a few hundred MeV, a range of values that not only can hardly be accepted, but 
is also such as to render the overall pionic corrections to the RFG responses quite 
negligible. It remains however to be verified whether the above findings are still 
valid in the 2p-2h sector or when mesons heavier than the pion are brought into 
play. 

As a general remark we note that, even in the presence of the vertex function 
F^, it is still possible to drop the 5 (smeared) piece of the OPEP while preserving 
gauge invariance via a suitable modification of the contact MEG current. Whether 
this should be done or not, however, is a question (as is the question about which 
value to choose for A^r) which should ultimately be answered by obtaining the best 
possible agreement with experiment. One might argue, however, that other mesons 
beyond the pion should be included in the formalism before attempting such a task: 
we intend to carry out this project in future work. 

In closing this section, we wish to address the question of the fc^-dcpendence of 
the various pionic contributions to the nuclear responses. In our view, kp is another 
parameter whose value should be set by comparing the theoretical predictions with 
experimental data. In figs. 17-20 we display the A;i?-dependence of the various 



contributions to the responses at q = 500 MeV/c divided by the functions 
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i^L(T,Vr;r7F) = .1 .: 3 ^ [(1 + ^)W2{r) - W,{t) + W2{t)A] (5.6a) 



and 



Krir, ijjr; Vf) = /^^"'^^ [2W^2(t) + W^2(t)A] (5.66) 

for the longitudinal and transverse channels, respectively. This permits us to focus 
only on the many-body part of the responses, and these reduced quantities will 
eventually scale for sufficiently high values of q. From fig. 17, where we display the 
self-energy contribution in the longitudinal electromagnetic channel, it is apparent 
that as kp increases this contribution has a different behaviour at high-a; than it 
does at low-w {i.e., for scaling variable ipr > and < 0, respectively — see ref. ^^) 
for a discussion of scaling in the RFG model). The net contribution to the total 
response, however, remains quite small. 

More pronounced is the sensitivity of the exchange terms to variations of kp, 
as illustrated in figs. 18a and b where the F and B longitudinal contributions are 
displayed. Indeed one sees that the F-term grows with hp and for l^/^^l ~ 0.8 
yields, in the longitudinal channel, a contribution as large as about 25% of the free 
response when kp = 250 MeV/c (corresponding nuclei around ^^Ca). The S-term 
also grows with kp, although its impact on the nuclear response reaches only about 
15% of the free response. As previously found, the transverse channel appears to 
be somewhat less affected by the F-correlations (fig. 19). Also the S-correlations 
are more significant in the longitudinal channel than in the transverse one and, 
overall, they are felt less than the F-correlations. The above results are clearly in 
line with the dominance of the (smeared) (5-term in shaping the nuclear response 
as previously discussed. 

As far as the MEG contributions (fig. 20) are concerned, their importance 
increases with hp more rapidly in the transverse channel than in the longitudinal 
one. For small kp {i.e., either light nuclei or regions of low density) they affect the 
nuclear responses only modestly; in particular, for = 225 MeV/c they contribute 
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only about 8% and 3% in the transverse and longitudinal channels, respectively. At 
higher kp (250 MeV/c) they may yield a contribution larger than 10%. 

As a last point to be discussed in this section, in parallel to the issue related to 
the impact of Att on the continuity equation discussed above, we briefly address the 
question of how global electromagnetic gauge invariance is affected by the size olkp- 
In this connection we consider the longitudinal/charge* nuclear response, since it is 
gauge invariance (through the continuity equation) that allows the replacement of 
the longitudinal component of the spatial electromagnetic current matrix elements 
with the charge matrix elements, or vice versa. We are not in a position to address 
this question at a quantitative level, because in the present paper, for sake of sim- 
plicity, we have ignored the convective current of the nucleon (only the dominant 
transverse spin current has been kept). We can, however, offer a simple guiding 
principle to help in recognizing whether gauge invariance is respected or not at the 
level of the many-body longitudinal/charge response: gauge invariance is fulfilled 
by a given Goldstone diagram contributing to the polarization propagator if its two 
electromagnetic vertices entail the same momentum flow. To understand what this 
means it helps to look at fig. 21. From the examples displayed there it is in fact 
apparent that, of the contributions calculated in the present paper, the self-energy 
diagrams respect gauge invariance, whereas the exchange ones do not. Analogously 
the antisymmetrized RPA (and even the ring diagrams alone. Fig. 21b) and ladder 
diagrams do not respect gauge invariance. 

These differences may be understood by considering the general expression for 

* The nomenclature used in studies of electroweak interactions with nuclei is not very 
good in this regard: what is referred to as the charge response in some places and the 
longitudinal response in others in fact usually involves both the = (charge) and fj, = 3 
or z (longitudinal) projections of the 4-current. As we discuss in this section, the two 
aspects are related by the continuity equation when the current is conserved and thus it 
has become common practice to use the terms interchangeably. Throughout this article 
we have used the word "longitudinal" to refer to all = and 3 aspects of electroweak 
current matrix elements. 
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the response written as the sum of three terms, namely 

-^j RL{q,oj) = R{QQ;q,uj) - —R,{QJ;q,uj) + -^R,,{JJ;q,uj) . (5.7) 

In the above R {QQ; q, u)), R^ {QJ', q, and R^z {J J', q, refer to responses with 
two charge, one charge and one longitudinal current and two longitudinal current 
electromagnetic vertices, respectively. The 2;-component of the nucleonic current, in 
the leading order of the non-relativistic reduction, involves the convection current 

Jz (x{2k + q)-q (5.8) 

and (2.1a), (3.1a) and (3.1b) contain integrals over k (actually, in the last two cases 
the integration is over a vector called ^2). Therefore, by exploiting the energy- 
conserving 5-function in the self-energy contribution (2.1a), it turns out that the 
current (5.8) can be expressed in this instance solely in terms of the external vari- 
ables q and u. Clearly, the same does not occur in the exchange contributions 
(3.1a) and (3.1b). Thus, via the continuity equation, for a translationally-invariant 
system treated at the level of the Hartree-Fock approximation, one has 

Rz{QJ;q,u;) = '^R{QQ;q,io) (5.9a) 



and 



^2 

Rzz{JJ;q,io)^ -^R{QQ-q,io). (5.96) 



Hence 

RLiq,uj)^R{QQ;q,uj) , (5.10) 

that is, RL{q,(^) may be expressed entirely in terms of a single response function. 
This is not the case, for example, in the case of the antisymmetrized RPA. It follows 
that for large kp, where the exchange contribution becomes relatively more impor- 
tant, it might even be necessary to account for the whole antisymmetrized RPA 
series, not to mention the ladder diagrams; in this instance serious consideration 
should be given to the gauge invariance of the theory. 
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6. Conclusions 

Several motivations prompted the present study: first, we wished to investigate the 
inclusive charge and spin electromagnetic nuclear responses within a model that as 
much as possible respects Lorentz and gauge invariance building on the work in 
ref. ^^). Secondly, our intention has been to explore, in parallel with the above, 
the weak neutral current responses (except for the axial-vector one whose analysis 
has been deferred to ref. ^^)). We have thus paved the way to the calculation of 
the asymmetry measured in parity-violating polarized electron scattering, which 
shall also be dealt with in ref. ^^). Finally our purpose has been to investigate the 
degree to which a specific nuclear model based on nucleonic and mesonic degrees 
of freedom, when treated consistently at the level of the forces and currents, yields 
meaningful physical results over a relatively large range of momentum transfers 
even when the framework assumed is limited to first-order perturbation theory. 

In this paper we have confined our attention to pionic effects, since the scale 
governing the range over which the currents and interactions are effective is espe- 
cially large for the pion, being the lightest meson. Thus we have started with pionic 
effects as our initial focus; this approach serves as a paradigm for the treatment of 
heavier mesons as well. Indeed, our basic approach in which it has been possible to 
explore (albeit, within a specific model of the nucleus) the interplay of forces and 
currents has natural extensions where heavier mesons can be included as well. 

Let us summarize our findings, beginning with the third of the motivations 
mentioned above. We have strived for consistency at two levels. The first one relates 
to the Feynman diagrams: of these we have retained all cases having one pionic line. 
This might appear insufficient for an adequate treatment of the forces and yet we 
have been able to show that results obtained in first-order perturbation theory are 
practically equivalent to those obtained in HF and RPA (especially to the first) when 
the force is carried by the pion. Accordingly, we are led to conjecture that even more 
sophisticated approaches, e.g., RPA with HF-dressed fermionic propagators, would 
not be substantially different from those using first-order perturbation theory (at 
least in the Ip-lh sector). This outcome is of importance because it would clearly 



be more diflficult to account for pionic MEC diagrams with two or more pionic lines. 

The second level of consistency relates to the continuity equation. Opera- 
tionally this involves asking whether the 5-function part of the OPEP, given its 
crucial impact on the nuclear responses, affects the continuity equation or not. 
The answer is no for pointlike nucleons and pions, and yes for extended ones. Ac- 
cordingly in the pointlike case one can, without affecting gauge invariance, drop 
the 5-function from the theory, arguing that short-range correlations will strongly 
quench its impact on the responses anyway. 

On the other hand, some aspects of this quenching might be phenomcnologically 
embedded in the factors Tt^. But then the resulting smeared 5 of the OPEP does 
contribute to the continuity equation. In such a case, to fulfill the latter one needs 
to modify the currents beyond the standard multiplicative factor F^, indeed by 
adding to the pion-in-flight current two extra pieces (that vanish when the mass 
term A^r oo). Their contribution to the nuclear responses in the Ip-lh sector is, 
however, quite small, although it should be checked whether this remains true in 
the 2p-2h sector and for heavier mesons. A final option, not explored in the present 
paper, namely to drop the smeared 5 of the OPEP, might deserve further study. 
In this case the restoration of gauge invariance can still be achieved by suitably 
modifying the contact current, namely, the only one contributed to by the 6. 

Turning to the second motivation, we have found a striking correlation effect 
in the longitudinal weak neutral current response R^y arising partly because the 
isovector and isoscalar components enter in R^y with a particular sign combination 
(dictated by the standard model), thus yielding a delicate cancellation in the case of 
the free RFG, and partly because the t — pionic correlations are much stronger, 
and of opposite sign, than the r = 1 ones. This outcome has an important bear- 
ing on the parity-violating polarized electron scattering asymmetry as discussed in 
ref. (see also ref. ^^)). 

Finally, coming to the first motivation above, as far as the forces are concerned 
we have found that: 

i) the charge and the spin responses are both hardened (shifted to higher a;) with 
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respect to the RFG; 

ii) the hardening is greater in the longitudinal/ charge channel (L) than in the spin 
channel (T); 

iii) the hardening fades away with increasing momentum transfer (at ~ 1 GeV/c 

very little is left of it) ; 
iiii) the hardening is due to the F^-modified 5-piece of the OPEP. 

Regarding the currents (MEG) we have found that: 

i) for the nuclear responses in the Ip-lh sector the MEG contributions are small, 
but not negligible (about 10% of the nuclear response in the transverse channel 
and at most 5% in the longitudinal one); 

ii) their contribution to the sum rules, relative to the RFG contribution, is almost 
constant as a function of q; 

iii) they always enter in such a way as to decrease the nuclear responses. 

In conclusion, our investigation of pionic effects started in ref. and pur- 
sued in the present paper shows the importance of the role they play in determin- 
ing the electromagnetic and weak neutral current quasielastic nuclear responses. 
Predictions regarding the g-dependent hardening of the longitudinal and trans- 
verse responses, Rl and Rt, as well as the dramatic modification of the longitudi- 
nal parity-violating response R^y emerge from our results. The approach taken, 
namely through extensions to the relativistic Fermi gas model where Lorentz co- 
variance and consistency in treating forces and currents can be maintained, allows 
us to explore the origin of these predictions in a relatively direct way. In partic- 
ular, it has proven to be instructive to see how they arise from the interplay of 
exchange and self-energy contributions, of central and tensor pieces of the force 
and of isoscalar/isovector correlation effects. Naturally, more remains to be done 
before a complete (hadronic) picture will be attained — in future work we intend 
to explore the nature of the quasielastic responses when heavier mesons and sum- 
mation schemes other than HF and RPA are incorporated using the same basic 
RFG-motivated framework. 
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Appendix A 

In this Appendix we perform the spin and isospin traces entering into the self-energy 
and exchange diagrams of fig. 4. 

a) Spin 

The spin matrix element to be calculated in the first-order self-energy term is 

Sl^T = 1] < Sa\Ol^T\Sb >< SbSc\K\SaSd >< Sd\dl^T\Sa > , {A.l) 

where 0£ j, is the spin operator associated to the 7NN vertex in the longitudinal 
or transverse channel, namely 

01 = 1 (^-2) 

and 

d^ = {ax g) . , (A.3) 

where i = 1,2 and 1 is the 2x2 unit matrix. In (A.l) can be either the central 
• (T2) or the tensor operator 

'S'i2(A;) = {Skikj - Sij) cruo-2j , {A A) 

k being the momentum carried by the pion. In the longitudinal channel one therefore 
gets 

St^central = X] < Sa\l\Sb >< Sb\ai\Sc >< Sc\cri\Sd >< Sd\l\Sa > 

= Tr{aiai} = 6, (A5) 

(repeated indices are meant to be summed), whereas the tensor force does not 
contribute to the longitudinal self-energy: 

SL%nsor = ^ < Sa\l\Sb >< Sb\cri\Sc >< Sc\crj\sd >< Sd\t\Sa > {Skikj - Sij) 

= Tr {(Tiaj} {Skikj - Sij) = 0. (^.6) 



In the transverse channel the central part of the pion-exchange potential yields 

^T^ntral = ^ < ^a\i^ ^ Q)k\sb >< Sb\cri\Sc >< Sc\cri\sd >< Sd\{(T X q)k\Sa > 

= Skim^krsQmQsTr {aiaiaiGr] = 12. (A. 7) 

Note that, in accord with the (A. 3) for Oj., although the sum over k in (A. 7) should 

be restricted to /c = 1, 2, it is unnecessary to take all of the components, since 

{(Txq)-{(Txq)^{ax q)i{a x g)i + (cr x q)2{(T x q)2 {A.8) 

if q points along the 2;-axis. Using similar arguments one can show that the tensor 
interaction gives no contribution to the self-energy in the transverse channel: 

qs.e. 

'^T,tensor 

< Salic X q)k\sb >< Sb\(7i\Sc >< Sc\crj\sd >< Sdliff X q)k\Sa > {Skikj - 6ij) 
= SklmSkrsqmqsi^kkj - Sij)Tr {(Ji(Ti(J jGr} = 0, (^.9) 

having used the property 

Tr {aiaiajar} = 2 (SuSjr + SirSij - SijSir) ■ (^.10) 
Let us now consider the exchange term. The spin factor is then given by: 



-lexch 



< Sa\Ol^T\Sb>< SbSd\Va\ScSa>< Sc\dl^j^\Sd> . {A.U) 

In the longitudinal response this leads to 

Sr,central = ^ < Sa\t\sb >< Sb\cri\Sc >< Sd\cri\Sa >< Sc\l\Sd > 

= Tr{aiai} = 6 (A12) 



and 



S 



T,tensor = XI < Sa\l\sb >< Sb\cri\Sc >< Sd\crj\Sa >< Sc\l\sd > {Skikj - 6ij) 



Tr {(Tiaj} (Skikj - Sij) = 0. 



(A13) 



In the transverse response 



cexch 
'~'T,central 



< Sa\{(T X q)k\Sh >< Sb\cri\Sc >< Sd\cri\Sa >< Sc\{<T X q)k\Sd > 



and 



= Sklm^krsqmQsTr {aiaiarCTi} = -4 



qexch 
T,tensor 



{A.U) 



< Sa\{(T X q)k\sb >< Sb\cri\Sc >< Sd\o-j\Sa >< Sc\{(T X q)k\sd > {Skikj - Sij) 
= ^klm^krsQmqaiSkikj - Sij)Tr {(J i(J i(T r(T j] =4 1 - 2, ■ qj . (A.15) 

Note that the transverse channel is the only one where the tensor interaction gives 
a nonzero contribution, namely, via the exchange diagram. 

b) Isospin 

Let us now consider the isospin matrix elements involved in the self-energy 
and exchange terms. At the 7NN vertex we have both isoscalar and isovector [z- 
component) dependences: for example, the charge and convection current operators 
involve 



= 

and the magnetization operator involves 

-1 + T, 



1 + 



{AAQa) 



O = 



2 
1 



{A.m) 



Here the 1-terms give rise to the isoscalar (r = 0) particle-hole responses, while 
the Tz-terms correspond to the isovector (r = 1) ones. 
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The self-energy diagram takes equal contributions from the traces in the two 
isospin channels: 

T'-^ {t=0) = j J2 < ta\t\tb >< tbtc\Ti ■ T2\tctd >< td\l\ta > 



4 

tbtctd 



= lTr{nn} = I {A.17) 



and 



■('^=1) = T XI < ta\rz\tb >< tbtcln ■ T2\tctd >< td\Tz\ta > 



rpS.e 

4 

tbtctd 

= \Tr{T,TiTiT,} = ^. (A.18) 

In contrast, in the exchange term the isoscalar and the isovector traces are different 
in magnitude and sign. They read 

T^^-'^(t = 0) = ^ J] < ta\t\tb >< tbtd\Ti ■ T2\tcta >< tc\l\td > 



4 



= \Tr{nn} = \ (^.19) 



and 

j^exch^^ = 1) = 1 < ta\rz\tb >< tbtd\ri ■ T2\tcta >< tc\T,\td > 

= ]TT{T,nT,n} = ~. (^.20) 

Finally let us consider the n^'^-order exchange diagram, which contains n pionic 
lines inside one bubble as in fig. 7. The corresponding isospin factor reads 

r^exch ^ Tr{d^TaTi2...T,n-inn0^nnnn-l-n2ni] {A.21) 

and can be obtained from the first-order ones by recursion relations: 

on 

T^-^'^(t = 0) = — (A.22) 



and 

C— 1 V 

T-^'^(r = 1) = -T-f (r = 1) = (A.23) 



Appendix B 

In this Appendix we first give the functions N'l,t and Tir entering in the formulae 
(3.1-4) for the exchange contributions to the response. The first two functions are 
the foUowing: 



X 



Att- 



ln|^ {ki,k2;z,ilj)\-ln\'iij-z\ 



+ 1 



(A^ + k! + klf - AklkliXl + kl + kl) - Ajj^kUXl -kl + kl) - 4i;zkKXl + kf - kj) 

[{XI + kl + klf - Aklkl]T~^ {kl, k2; z, ^) 

+ {kiM; z, V) - m\- {ki,k2; z, i;) (5.1) 

and 

k2; z, ^) = I ^\n\2klz - ^Xl + /c? + kl) + kiJ^~^ {ku k2; z, 

^ {XI + kl + klf - A{Xl + kl + kl){klkl + Xl^z) + Atl;kl{tl; - z){Xl - kj + kl) 
[{XI + kj + klf - Aklkl]T~ {kl, k2; z, ^) 

+ Mp- {kuk2; z, i;) - M^^ {ki,k2; z, V^), {B.2) 

where 

Ml{kuk2; z, iP) = - ^^f^ [ln|6?A(A:i, k2; z,^)\- - z\] 

1 k2 

MT{ki,k2;z,'4)) = -^ — [Tx{ki,k2;z,'4)) 

+-^(A2 -kl + kl)\n\2klz - V'(A2 + kl + kl) + kiJ^x{ki, k2; z,il^)\] 
kl 

Q\{ki,k2; z,il)) = 2i(;{i(; - z){X'^ - kl + kl) + Tx{ki, k2; z,il))T\{ki,k2;ii,il^) 
+ Tl{ki,k2;'il),ip) 



Tx{ki,k2; z, i;) = {X^ + kl + kl)^ - #2(A2 + /cf + kl) + Aklz^ + 4/^1^^ - ^klk 

(S.3) 

We have also set Att = A^^/kp and = rriT^/kF. The function Ht, which gives in 
a fully analytic form the transverse exchange response at q > 2kF, reads 
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Xl + {k-p+lf 



Xl + {k + p-l)^ 



{BA) 



where 



T\k,p) = 



1 



1 



tan 



1 fk -p+l 
A 



tan 



^ [k + p-l 



X 



+ - 6/c2(l - A^) - 3(1 + A2)2 + 12p(A2 + fc2 + 1) 

A^ +(^k-p + lf 



6p^(A^ + /c^ + 3) - 3p^(p - 4)]ln 



A2 + (/c + J9 - 1)2 



+- (p - l)3ln [A^ + (/c - p + 1)2] [A^ + (fc + p - 1)^1 



(S.5) 



Finally, wc give the expression for the real part of the first-order forward-going 
(F) longitudinal polarization propagator (see eq. (3.6)), limited for simplicity to the 
case q > 2k p'- 



1 



ml + {ki - 



1 



Lo — \Q'^\/2'mN — q ■ ki/niN ^ — IQ^I/Smjv — 9 ■ k2/'mN 

- 5{uj - \Q^\/2mN - q ■ fci/wjv) - \Q'^\/2mN - q ■ k2/mN) 



{B.6a) 



X 



1287r^ r/Fn'^ rriN 

/ dki / dxi / dk2[V{ki,k2;kiXi,'i(;r) -V{ki,-k2;kixi,'i(;r)] 
Jo J-i Jo 



+ 77^/ dkkln 



1 + ml-k^ + -ml- fc2)2 + 4m2 (1 - V'r) 



1 + A2 - A;2 + ^ (1 - A2 - A;2)2 + 4A2 (1 - 



XI - ml TT^ 



--{Xl + 2- 2i^l - A.^A2 + 4 - 4V^2) 

2 



{EM) 
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where Xr=Q = 3/2, Ir=i = -1/2, 



X {inlV' + k2\ + ln|2yi(7A - k2){X^ - kf + k^) - AkHiJ: - k^^iP - yi) 
+ (A^ + kl + kl- 2yik2)J^x{ki, k2; yi,i^) + J^l{ki, /ca; yi» | }• 

{B.7b) 
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Figure Captions 
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Fig. 1. Goldstone diagrams for the self-energy, with the pion dressing the 
particle (a) or the hole (b) line. 

Fig. 2. The self-energy contribution to the longitudinal response function is 
shown as a function of the energy transfer u for three values of the momen- 
tum transfer: g=300, 500, and 1000 MeV/c. The dashed lines represent the 
separated particle and hole terms, the solid line being their sum. Here and in 
the following figures, unless otherwise specified, the nucleus considered is ^^C, 
/CF=225 MeV/c and A^=1.3 GeV. 

Fig. 3. The Hartree-Fock longitudinal response (solid curve) is displayed 
and compared to the free RFG (dotted) and first-order self-energy correlated 
(dashed) responses. The first and last almost coincide. 

Fig. 4. Goldstone diagrams corresponding to the forward-going (a) and backward- 
going (b) terms of the exchange correlations, defined in eqs. (3.3) and (3.4). In 
the text these are labeled "F" and "S" , respectively. 

Fig. 5. The exchange contribution to the longitudinal response is shown (solid 

curves). The forward- (dashed) and backward-going (dot-dashed) terms, "F" 
and "S" , respectively, are also separately displayed. 

Fig. 6. The same as fig. 5, but now for the transverse response. 

Fig. 7. Goldstone diagrams entering into the Tamm-Dancoff series for the 
polarization propagator with pion exchange. 

Fig. 8. The free (dotted), first-order exchange (dashed) and first continued 
fraction (solid) longitudinal responses are plotted as functions of u at ^=500 
and 1000 MeV/c. For the results in the left-hand panel, the interaction is the 
pure "(5-like" part of the pion potential, while for the results in the right-hand 
panel the full potential is taken. 
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Fig. 9. The free Fermi gas longitudinal response is shown as dotted curves; the 
dashed curves include the self-energy and exchange correlations, whereas the 
solid curves also include MEC. 

Fig. 10. The same as fig. 9, but now for the transverse response. 

Fig. 11. The shift of the quasielastic peak due to pionic correlations is shown as 
a function of q in the longitudinal (solid curve) and transverse (dashed curve) 
channels. 

Fig. 12. The longitudinal (upper panel) and transverse (lower panel) sum rules 
Sl{(i) and St{q), as defined through eq. (4.5), are shown as functions of q. 
The dotted curves correspond to the RFG, the dashed curves include self- 
energy and exchange correlations and the solid ones contain as well the MEC 
contribution. In the small windows the ratios of the correlated sum rules to 
the free ones are displayed: the dashed curves do not include the MEC, the 
solid ones do. 

Fig. 13. Upper panel: the isoscalar free (dashed) and correlated (solid) lon- 
gitudinal responses are displayed as functions of uj at g=500 MeV/c. Lower 
panel: the same quantities are shown in the isovector channel. 

Fig. 14. The longitudinal weak neutral current response R^y at q= 500 MeV/c, 
showing the free RFG (dashed) and correlated (solid) results. 

Fig. 15. The exchange contribution to the longitudinal response at q=500 
MeV/c is displayed in the upper panel for three different values of A^r: 10 GeV 
(dashed), 1.3 GeV (dot-dashed) and 0.8 GeV (dotted). In the lower panel the 
longitudinal self-energy is shown, with the same meaning for the three curves. 

Fig. 16. The exchange (upper panel) and MEC (lower panel) contributions 
to the transverse response at q=500 MeV/c are displayed for At^= 10 GeV 
(dashed), 1.3 GeV (dot-dashed) and 0.8 GeV (dotted). 

Fig. 17. The self-energy contribution to the longitudinal response divided by 
the function Kl{t, ipr'i Vf) (eq. 5.6a) at q'=500 MeV/c is displayed as a function 
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of the scaling variable ipr', three different values of kp have been chosen: 200 
MeV/c (dashed), 225 MeV/c (dot-dashed) and 250 MeV/c (solid). 

Fig. 18. The forward-going (F, upper panel) and backward-going {B, lower 
panel) terms of the exchange, divided by KLir^'^riVp) in the longitudinal 
channel shown as functions of V'r for ^^=200 MeV/c (dashed), 225 MeV/c 
(dot-dashed) and 250 MeV/c (sohd) at q=500 MeV/c. 

Fig. 19. The forward-going (F, upper panel) and backward-going (S, lower 
panel) terms of the exchange, divided by Kt{t, i/jr', Vf) (eq. 5.6b) in the trans- 
verse channel are shown as functions of t/jr for /cf=200 MeV/c (dashed), 225 
MeV/c (dot-dashed) and 250 MeV/c (sohd) at g=500 MeV/c. 

Fig. 20. The MEC contribution to the longitudinal response divided by Kl{t, ipr', Vf) 
is displayed in the upper panel as a function of if^r for /cf=200 MeV/c (dashed), 
225 MeV/c (dot-dashed) and 250 MeV/c (solid) at q=500 MeV/c. In the lower 
panel the transverse MEC contribution divided by Kxij^'^'r'i'nF) is shown for 
the same set oi kp values. 

Fig. 21. Gauge invariance in many-body theories. Diagram (a) (self-energy) 
has the same momentum flow in the two electromagnetic vertices. Diagrams (b) 
(ring), (c) (exchange), (d) (exchange contributions to antisymmetrized RPA) 
and (e) (ladder) have not. The first class of diagrams fulfills gauge invariance, 
the second does not. 



